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ABSTRACT: SARS (severe acute respiratory syndrome) is caused by a newly discovered coronavirus. A 
key enzyme for the maturation of this virus and, therefore, a target for drug development is the main 
protease 3CL?™ (also termed SARS-CoV 3CL?"°). We have cloned and expressed in Escherichia coli the 
full-length SARS-CoV 3CL? as well as a truncated form containing only the catalytic domains. The 
recombinant proteins have been characterized enzymatically using a fluorescently labeled substrate; their 
structural stability in solution has been determined by differential scanning calorimetry, and novel inhibitors 
have been discovered. Expression of the catalytic region alone yields a protein with a reduced catalytic 
efficiency consistent with the proposed regulatory role of the a-helical domain. Differential scanning 
calorimetry indicates that the a-helical domain does not contribute to the structural stability of the catalytic 
domains. Analysis of the active site cavity reveals the presence of subsites that can be targeted with 
specific chemical functionalities. In particular, a cluster of serine residues (Ser139, Ser144, and Ser147) 
was identified near the active site cavity and was susceptible to being targeted by compounds containing 
boronic acid. This cluster is highly conserved in similar proteases from other coronaviruses, defining an 
attractive target for drug development. It was found that bifunctional aryl boronic acid compounds were 
particularly effective at inhibiting the protease, with inhibition constants as strong as 40 nM. Isothermal 
titration microcalorimetric experiments indicate that these inhibitors bind reversibly to 3CL?° in an 
enthalpically favorable fashion, implying that they establish strong interactions with the protease molecule, 


thus defining attractive molecular scaffolds for further optimization. 


From November 2002 through June 2003, more than 8000 
persons contracted severe acute respiratory syndrome (SARS), 
a highly contagious respiratory infection never seen before 
in humans. The overall mortality rate associated with SARS 
has been estimated to be close to 10% by the World Health 
Organization, but exceeded 50% in patients at least 65 years 
of age. In March 2003, 4 months after the first case was 
reported, a new coronavirus was identified as the etiological 
agent of SARS. In April 2003, the complete genome of the 
virus was sequenced, providing the opportunity for more in- 
depth studies and the identification of potential targets for 
drug and vaccine development. The genome sequence of the 
virus immediately revealed the presence of a protease (3CL?” 
or SARS-CoV 3CL?°) whose sequence was highly identical 
to those of proteases known to be essential for the life cycle 
of other coronaviruses (refs /—3 and unpublished data from 
this laboratory). In July 2003, one crystallographic structure 
of 3CL?° was released (PDB entry 1q2w by Bonnano et al.). 
A second structure was published in November 2003 (4). 


The active site of the SARS-CoV 3CL? contains a 
catalytic dyad defined by His41 and Cys145, which is similar 
to the arrangement found in other coronavirus proteases. 
From a structural point of view, the catalytic dyad is located 
in a predominantly /-structure (residues 1—184) reminiscent 
of the two-domain fold found in the chymotrypsin family 
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FiGurRE 1: Structure of 3CL?"° (PDB entry 1q2w) showing the 
catalytic domains in red and the a-helical domain in blue. Catalytic 
residues His41 and Cys145 are shown using a van der Waals 
representation. The color transition (Gln192-Ala193) indicates the 
point at which a stop codon was placed to express the truncated 
version of 3CL?, 


(Figure 1). Residues 201—303 form a third compact a-helical 
domain connected to the catalytic 6-domains by a long loop 
region (residues 185—200). The precise biological role of 
the a-helical domain is not completely understood, even 
though it has been hypothesized that it might play a role in 
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substrate recognition in other coronaviruses (5, 6). From the 
dimeric structure of SARS-CoV 3CL?”, Yang et al. (4) have 
concluded that the helical domain might be responsible for 
orienting the amino terminus of one monomer to a position 
in which it can interact with the active site of the other 
monomer in the dimer. In solution, SARS-CoV 3CL?" is 
characterized by a weak dimerization constant of 100 uM 
which is equivalent to 3.4 mg/mL (7). 

For the studies presented here, we have cloned and 
expressed in Escherichia coli the full length SARS-CoV 
3CL?" as well as a truncated form containing only the 
catalytic domains. The recombinant proteins have been 
characterized enzymatically using fluorescently labeled 
substrates; their structural stability in solution has been 
determined by differential scanning calorimetry, and novel 
inhibitors have been discovered. Further characterization of 
these inhibitors suggests that they provide good scaffolds 
for further optimization and drug development against 
SARS. 


EXPERIMENTAL PROCEDURES 


Cloning. cDNA corresponding to the SARS 3CL?® gene 
(Tor2 strain, GenBank entry AY274119) inserted into a 
pBR194c vector was kindly provided by the British Colum- 
bia Cancer Agency Branch (Vancouver, BC). Competent 
cells (XL-1 Blue, Stratagene) were transformed for plasmid 
propagation under ampicillin selection. Plasmid DNA was 
isolated (Plasmid Midi Kit, Qiagen), and the gene was 
amplified by PCR with appropriate primers using Pfu Turbo 
DNA Polymerase (Stratagene). To prevent artifacts, the 
original plasmid was degraded with a DpnlI (Stratagene) 
digestion reaction (1 h at 37 °C) followed by inactivation of 
DpnI (20 min at 80 °C). Blunt-end directional cloning was 
performed by a topoisomerase reaction into a pET100 vector 
(Champion pET Directional TOPO Expression and Cloning 
Kit, Invitrogen). The protease gene was cloned in frame with 
an N-terminal peptide containing a polyhistidine tag for 
further purification by affinity chromatography. An enter- 
okinase recognition site was present for removal of the 
amino-terminal tag after purification. As a result, the 
recombinant protease contains five additional residues (Asp, 
His, Pro, Phe, and Thr) at the amino terminus. One Shot 
TOP10 competent cells (Invitrogen) were transformed with 
the product from the cloning reaction. Plasmid DNA was 
purified, and the gene insertion and its directionality, as well 
as the integrity of the pET100 vector, were confirmed by 
DNA sequencing. BL21 Star DE3 (Invitrogen) competent 
cells were transformed for protein expression under ampi- 
cillin selection and IPTG induction. 

The gene encoding the truncated protein containing only 
the catalytic domains was obtained by placing a stop codon 
in the wild-type plasmid at the position corresponding to 
Alal93 (Quick Change Site Directed Mutagenesis Kit, 
Stratagene). After transformation as described above, the 
plasmid DNA was isolated (Plasmid Midi Kit, Qiagen) and 
the mutation confirmed by sequencing. BL21 Star DE3 
(Invitrogene) competent cells were transformed for protein 
expression under ampicillin selection and IPTG induction. 
The truncated protein was expressed and purified in the same 
manner as wild-type SARS 3CL?”. 

Protein Expression and Purification. Plasmid-encoded 
SARS 3CL? was expressed as a soluble fraction in BL21 
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Star DE3 E. coli competent cells (Invitrogen). Cells were 
grown in LB supplemented with ampicillin (50 “g/mL) at 
37 °C, induced with IPTG when the optical density was 
~0.8, and harvested after 4 h. Cells were resuspended in 
lysis buffer [50 mM potassium phosphate (pH 7.8), 400 mM 
sodium chloride, 100 mM potassium chloride, 10% glycerol, 
0.5% Triton-X, and 10 mM imidazole] and broken with three 
passes through a French pressure cell (216000 psi). Cell 
debris was collected by centrifugation (20000g at 4 °C for 
20 min). The supernatant was filtered using a 0.22 um pore 
size filter (Millipore) and applied directly to a nickel affinity 
column (HisSelect, Sigma) that had been pre-equilibrated 
with binding buffer [50 mM sodium phosphate, 0.3 M 
sodium chloride, and 10 mM imidazole (pH 8.0)]. The 
protease was eluted with a linear gradient of elution buffer 
[50 mM sodium phosphate, 0.3 M sodium chloride, and 250 
mM imidazole (pH 8.0)]. During the concentration of 
protease fractions, the elution buffer was exchanged gradually 
with storage buffer [10 mM sodium phosphate, 10 mM 
sodium chloride, 1 mM TCEP, and 0.5 mM EDTA (pH 7.4)]. 
The polyhistidine tag of the fusion protein was cut through 
incubation with 0.1 unit of enterokinase (Invitrogen) per 56 
ug of protease for 48 h at 4 °C. The required number of 
units was determined by setting up pilot reactions using 
different amounts of enterokinase at different temperatures, 
and then verifying the efficiency of the reaction by SDS— 
PAGE. This reaction mixture was passed again through the 
nickel affinity column, and the flow-through containing the 
protease was collected. Any residual enterokinase was 
removed by incubation with EK-away resin (Invitrogen). The 
sample was diluted 1:4 with storage buffer and concentrated 
(> 10 mg/mL). The purified protein was then stored at —20 
°C. The sample was more than 95% pure, as assessed by 
SDS—PAGE. 


Enzymatic Activity. The activity of the full-length and 
truncated 3CL?"? was measured by continuous kinetic assays 
using the commercially available fluorogenic substrate 
Dabcyl-Leu-Ala-Gln-Ala-Val-Arg-Ser-Ser-Ser-Arg-Edans 
(Bachem). The fluorescence intensity was monitored in a 
Cary Eclipse fluorescence spectrophotometer (Varian) using 
wavelengths of 330 and 500 nm for excitation and emission, 
respectively. The experiments were performed with the same 
buffer used to store the enzyme [10 mM sodium phosphate, 
10 mM sodium chloride, 1 mM TCEP, and 0.5 mM EDTA 
(pH 7.4)]. Kinetic parameters, K, and Kat, were determined 
by initial rate measurements at 25 °C. The reaction was 
initiated by adding protease (final concentration of 1 4M) 
to a solution containing different final concentrations of the 
substrate (O—100 uM). 


Inhibition Kinetics. FL compounds for inhibition studies 
were kindly supplied by Fulcrum Pharmaceuticals (Larch- 
mont, NY). Inhibition constants and the inhibition mechanism 
were assessed by determining the apparent kinetic parameters 
at (1) a constant substrate concentration and different 
inhibitor concentrations and (2) constant inhibitor concentra- 
tions and variable substrate concentrations. Protease (final 
concentration of | “~M) was incubated for 10 min at room 
temperature with the inhibitor (final concentration of 0—180 
UM) and the reaction initiated by adding substrate to the 
desired final concentration (5 uM for experiments at a 
constant substrate concentration). 
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The dependence of the enzymatic rate in terms of substrate 
and inhibitor concentrations was analyzed in terms of the 
general equation (see ref 8) 


vy K,, + [IS] 


Uo Kt 4 St a 


1 1 


where vy and vo are the initial rates at a given inhibitor 
concentration and without inhibitor, respectively, [S] is the 
concentration of the substrate, [I] is the concentration of the 
free inhibitor, K,, is the Michaelis constant, and K; is the 
inhibition constant. The parameter a equals co for competitive 
inhibition, approaches zero for uncompetitive inhibition, and 
assumes a finite value for noncompetitive inhibition (8). All 
the parameters (Ky, Ki, and &) were determined by global 
nonlinear regression analysis. 

Isothermal Titration Calorimetry. Isothermal titration 
calorimetry experiments were carried out using a high- 
precision VP-ITC titration calorimetric system (Microcal 
Inc.). The enzyme solution in the calorimetric cell was titrated 
with inhibitor solutions dissolved in the same buffer [10 mM 
sodium phosphate, 10 mM sodium chloride, | mM TCEP, 
and 0.5 mM EDTA (pH 7.4)]. The heat evolved after each 
ligand injection was obtained from the integral of the 
calorimetric signal. The heat due to the binding reaction 
between the inhibitor and the enzyme was obtained as the 
difference between the heat of reaction and the corresponding 
heat of dilution. 

Differential Scanning Calorimetry. The heat capacity of 
full-length and truncated 3CL?"® was measured as a function 
of temperature with a high-precision differential scanning 
VP-DSC microcalorimeter (Microcal Inc., Northampton, 
MA). Protein samples and reference solutions were properly 
degassed and carefully loaded into the cells to avoid bubble 
formation. Thermal denaturation scans were performed with 
freshly prepared buffer-exchanged protease solutions in 10 
mM sodium phosphate, 10 mM sodium chloride, 1 mM 
TCEP, and 0.5 mM EDTA (pH 7.4). Data were analyzed by 
software developed in this laboratory. 


RESULTS AND DISCUSSION 


Enzymatic Activity Measurements. A fluorescence-based 
enzyme activity assay was implemented to evaluate the 
inhibitory potency of potential inhibitors. A commercially 
available fluorogenic peptide (Dabcyl-Leu-Ala-Gln-Ala-Val- 
Arg-Ser-Ser-Ser-Arg-Edans) that contains the SARS 3CL?° 
Gln-Ala cleavage sequence was used to measure the catalytic 
rate constants (Kea) and the Michaelis constants (Ky) as 
described in Experimental Procedures. Figure 2 shows the 
dependence of the initial enzyme rate on the fluorogenic 
substrate concentration. The enzyme velocity data are 
consistent with a Km of 9 + 1 uM and a catalytic constant 
(Kear) of 0.033 + 0.001 s~ for full-length 3CL"° at 25 °C. 
Incidentally, the N-terminal His-tagged protease exhibited 
similar enzymatic activity, indicating that the addition of the 
six histidyl residues to the amino terminus does not 
significantly affect catalysis. The Km of the fluorogenic 
substrate is significantly lower than those reported for 
unlabeled undecapeptides with sequences corresponding to 
the cleavage sites in the SARS coronavirus polyprotein (7). 
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FiGure 2: Enzymatic activity of full-length (@) and truncated (O) 
3CL?. The truncated 3CL? contains only the catalytic domains 
(residues 1—192). Shown in the figure is the initial velocity in 
arbitrary fluorescence units as a function of substrate concentration. 
The experiments were performed at 25 °C in 10 mM sodium 
phosphate, 10 mM sodium chloride, | mM TCEP, and 0.5 mM 
EDTA (pH 7.4). 


The fluorogenic substrate contains the correct Gln-Ala 
cleavage sequence but differs from the viral consensus 
sequences at positions P4 and P2 [P2, however, has a lower 
specificity (4)]. The lower K,, for the fluorogenic substrate 
is indicative of a higher affinity for the protease. It is not 
clear whether the K,, differences are real or the result of 
different experimental conditions, protein preparations (e.g., 
N-terminal extension or C-terminal His tag), or assay 
methods. In fact, similar differences have also been observed 
for spectroscopic and HPLC-based enzymatic assays of 
HIV-1 protease substrates (9—2/). 

Since the helical domain of 3CL? is supposed to play a 
role in substrate recognition (5, 6), we decided to investigate 
the enzyme activity of the isolated catalytic domains. Toward 
that end, a truncated protein containing only residues 1— 192 
(shown in red in Figure 1) was also expressed. The catalytic 
domains alone also exhibited enzymatic activity (Figure 2), 
albeit with a reduced catalytic efficiency. Its activity is 
characterized by a Ky of 13 + 2 uM and a Kea of 0.014 + 
0.001 s~!. The catalytic efficiency of the full-length protease 
(Kea/Km) is ~3.5 times higher than that of the truncated 
protein. Previously, Anand et al. (6) reported an even larger 
loss in activity for another coronavirus protease (TGEV 
3CL?*). These differences and the reported preference for 
substrates with /-sheet-like structures (7) suggest that, at a 
quantitative level, the regulatory role of the helical domain 
in correctly positioning the substrate for catalysis (4—6) 
might depend on the size and structural characteristics of 
the substrate. 

Structural Stability of SARS 3CL?". Figure 3 shows the 
temperature dependence of the heat capacity function of full- 
length 3CL?”. The thermal denaturation transition is centered 
at 52.5 °C under the conditions of these experiments. At pH 
7.4, the transition is followed by precipitation, leading to a 
lack of reversibility in the microcalorimetric scans that 
precludes a full thermodynamic analysis. Nevertheless, it can 
be concluded that the transition is highly cooperative and 
proceeds in a two-state fashion with no indication of 
independent behavior by either of the domains. Experiments 
performed at concentrations ranging between 0.17 and 0.4 
mg/mL did not exhibit an upward trend in denaturation 
temperature as could be expected for a dimeric protein. 
However, these concentrations are well below the dimeriza- 
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FiGURE 3: Temperature denaturation of full-length (A) and 
truncated (B) 3CLP° as determined by differential scanning 
calorimetry. Calorimetric scans for both proteins were performed 
at identical protein concentrations (0.4 mg/mL) at a scanning rate 
of 1 °C/min. 


tion Kg, and the precipitation observed after denaturation will 
mask any such dependence. Analytical gel filtration experi- 
ments by Fan et al. (7) indicate that the dimer dissociation 
midpoint occurs at 3.4 mg/mL. Analytical ultracentrifugation 
experiments (not shown) are consistent with those results 
and indicate that at pH 7.4 and protein concentrations of <1.5 
mg/mL the protease is mainly monomeric. Microcalorimetric 
experiments were also performed with the truncated 3CL””° 
containing only the catalytic domains (Figure 3). The 
denaturation transition of the truncated 3CL?” is centered at 
53.3 °C, a temperature slightly higher than that of full-length 
3CL?"?, and is also consistent with two-state behavior. In this 
respect, it has been known for many years that the structurally 
similar protease chymotrypsin obeys two-state behavior 
despite being composed of two structural domains similar 
to those of 3CL*° (22). The observation that the truncated 
protein has a slightly higher stability than the full-length 
protein is an indication that the a-helical domain is intrinsi- 
cally less stable and that it derives a significant fraction of 
its stabilization energy from its interactions with the catalytic 
domain (see, for example, refs 23 and 24). 

Inhibitor Design. Catalytic residues His41 and Cys145 are 
located in a relatively shallow cavity between the two 
domains defining the catalytic site (Figure 4). Close to the 
catalytic residues is a cluster of serine residues (Ser139, 
Ser144, and Ser147) that define an attractive subsite for the 
design of high-affinity inhibitors. Serl144 and Serl47 are 
located in a cavity that can be targeted by low-molecular 
weight compounds. The serine cluster is highly conserved 
in all known coronavirus proteases as illustrated in Figure 
5, suggesting that targeting this site might lead to the 
development of wide-spectrum coronavirus protease inhibi- 
tors. In particular, Ser139 and Ser147 are conserved in all 
known coronaviruses. Ser144 is conserved in 55% of the 
reported sequences, while in an additional 15% of the 
sequences, the closely related amino acid threonine occupies 
that position. The remaining 30% are occupied by alanine. 
Targeting a highly conserved region among different coro- 
naviruses provides two advantages, the potential for wide- 
spectrum antivirals and better inhibitor resilience against 
mutations associated with drug resistance. 

Because of the known potential reactivity of boronic acid 
compounds with the hydroxyl group in serine residues, the 
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FIGURE 4: Catalytic site of 3CL?° showing the location of Ser139, 
Serl44, and Ser147 (red). Catalytic residues His41 and Cys145 are 
shown in yellow. The serine cluster defines an attractive subsite 
for targeting by potential inhibitors. 


inhibitory potency of bifunctional boronic acid compounds 
was evaluated. A chemical scaffold containing two phenyl 
boronic groups (FL-078, chemical structure given in Table 
1) attached to a central aromatic ring by ester linkages was 
tested. The compound was shown to have inhibitory activity. 
Figure 6A shows the dependence of enzyme activity on the 
concentration of compound FL-078. The inhibitor is char- 
acterized by an inhibition constant in the low micromolar 
range (Kj = 4.5 uM). To elucidate the inhibition mechanism 
of this chemical scaffold, inhibition kinetic experiments were 
performed at different constant inhibitor concentrations. The 
results are shown in Figure 6B and indicate that the inhibitor 
affects not only the apparent K, but also the Vinax. This 
behavior is consistent with a molecule that binds to both the 
free enzyme and the enzyme—substrate complex. The 
reversibility of the inhibition mechanism was examined by 
incubating the enzyme with excess inhibitor (200 uM) and 
then diluting the inhibitor concentration immediately before 
the enzyme activity assay. The kinetic parameters were the 
same in both cases. Accordingly, the data were analyzed by 
global nonlinear least-squares analysis with the general 
equilibrium model discussed by Copeland (8) and the 
equations described in Experimental Procedures. Different 
variations of FL-078 were tested, including several isomers, 
replacement of the central aromatic ring with shorter linkers, 
different functionalities at the linker region, and different 
functionalities at the phenyl boronic rings. The highest 
improvement in affinity was observed when the ester linkage 
between aromatic rings was replaced with an amide group, 
resulting in nanomolar binding affinities. Table 1 summarizes 
the inhibition constants (Ki) for five representative com- 
pounds. Besides the inhibition constants, of particular interest 
is the parameter a, which ranges between 1.5 and 5.6 for 
the inhibitors included in Table 1. For a purely competitive 
inhibitor, © = co, whereas for a noncompetitive inhibitor, a 
ranges somewhere between | and . If a = 1, the inhibitor 
binds with equal affinity to the free enzyme and to the 
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FiGureE 5: Sequence alignment of the region containing the serine cluster (Ser139, Ser144, and Serl147) among proteases from different 
coronaviruses. GenBank accession numbers for protein sequences follow: NP_828863 for SARS-HCV (severe acute respiratory syndrome 
human coronavirus), AAF69341 for MHV ML-10 (murine hepatitis virus strain ML-10), NP_740610 for MHV A59 (murine hepatitis virus 
strain A59), P19751 for MHV JHM (murine hepatitis virus strain JHM), AAF19383 for MHV-2 (murine hepatitis virus strain 2), AAF69331 
for MHV Penn 97-1 (murine hepatitis virus strain Penn 97-1), AAF68919 for MHV ML-11 (murine hepatitis virus strain ML-11), AAL40396 
for BCV Quebec (bovine coronavirus strain Quebec), AAL57315 for BCV LUN (bovine coronavirus strain LUN), AAA64744 for BCV 
Mebus (bovine coronavirus strain Mebus), NP_742132 for BCV ENT (bovine coronavirus strain ENT), NP_839959 for PEDV-CV777 
(porcine epidemic diarrhea virus strain CV777), CAA83979 for TGEV Purdue-115 (transmissible gastroenteritis virus strain Purdue-115), 
NP_840003 for TGEV PUR46-MAD (transmissible gastroenteritis virus strain PUR46-MAD), AAK09095 for FIPV 79-1146 (feline infectious 
peritonitis virus strain 79-1146), NP_835346 for HCV 229E (human coronavirus strain 229E), NP_740623 for IBV Beaudette (avian infectious 
bronchitis virus strain Beaudette), AAQ21584 for IBV LX4 (avian infectious bronchitis virus strain LX4), CAC39112 for IBV Beaudette 
CK (avian infectious bronchitis virus strain Beaudette CK), and AAP92674 for IBV BJ (avian infectious bronchitis virus strain BJ). 
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FiGuRE 6: Inhibition of SARS-associated 3CL?° by compound FL-078 (see Table 1). Shown in panel A is the decrease in enzyme activity 
as a function of inhibitor concentration at 5 ~M S. Shown in panel B is the dependence of enzyme activity on substrate concentration at 
different constant inhibitor concentrations. Analysis of the data indicates that the inhibitors considered here are noncompetitive and affect 


both Ky, and Vinax (see the text for details). Solid lines correspond to a best parameter fit (Table 1). 


enzyme—substrate complex. For the inhibitors in Table 1, 
the affinity is higher for the free enzyme by factors ranging 
between 1.5 and 5.6, but significant affinity for the enzyme— 
substrate complex is maintained. This is somewhat expected 
from the structural location of the targeted serine cluster in 
relation to the catalytic dyad. It can be suggested, however, 
that inhibition of substrates larger than the decapeptide 
substrate used in these studies might exhibit a more competi- 
tive nature if they require a larger binding footprint and their 
binding is additionally compromised by the presence of the 
inhibitor. Analysis of the inhibition data in terms of Dixon 


plots (not shown) yielded a linear dependence of 1/v on 
inhibitor concentration (r? = 0.97), indicating that the 
enzyme—substrate—inhibitor complex is not catalytically 
active; i.e., the parameter / is equal to zero according to 
standard nomenclature (8, 25). The truncated form of the 
enzyme was equally inhibited by the bifunctional aryl boronic 
acid compounds studied here, consistent with the idea that 
these compounds target the catalytic and not the a-helical 
domain. 

Isothermal titration calorimetric experiments (Table 1) 
indicate that the inhibitors studied here bind to 3CL” with 
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Table 1: 3CL?® Inhibition Parameters for Bifunctional Aryl Boronic Compounds 


Ki AH -TAS 


favorable enthalpic and entropic contributions. This ther- 
modynamic profile is characteristic of a binding process 
driven by favorable inhibitor—protein interactions and not 
only by water exclusion due to hydrophobicity (26—28). It 
has been argued that enthalpically favorable lead compounds 
must be preferred over entropically driven ones since further 
optimization can be more readily accomplished (26-28). 

The compounds studied here have molecular weights 
ranging between 344 and 524 and are characterized by 
inhibition constants as low as the midnanomolar range; as 
such, they provide attractive scaffolds for drug design. 
Currently, we are engaged in the optimization of these 
scaffolds as potential anti-SARS drugs. The results will be 
the topic of a future communication. 
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